Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



o 
o 

o 

Q 

(N 



Fermi liquid theory for the nonequihbrium Kondo effect at low bias voltages 
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In this report, we describe a recent development in a Fermi liquid theory for the Kondo effect 
in quantum dots under a finite bias voltage V. Applying the microscopic theory of Yamada and 
Yosida to a nonequilibrium steady state, we derive the Ward identities for the Keldysh Green's 
function, and determine the low-energy behavior of the differential conductance dl /dV exactly 
up to terms of order {eV)'^ for the symmetric Anderson model. These results are deduced from 
the fact that the Green's function at the impurity site is a functional of a nonequilibrium dis- 
tribution /off((x)), which at eV = coincides with the Fermi function. Furthermore, we provide 
an alternative description of the low-energy properties using a renormalized perturbation the- 
ory (RPT). In the nonequilibrium state the unperturbed part of the RPT is determined by 
the renormalized free quasiparticles, the distribution function of which is given by /eff(a;). The 
residual interaction between the quasiparticles U , which is defined by the full vertex part at 
zero frequencies, is taken into account by an expansion in the power series of U . We also discuss 
the application of the RPT to a high-bias region beyond the Fermi-liquid regime. 

KEYWORDS: Kondo effect, Fermi liquid, Nonequilibrium, Keldysh formalism, Anderson model. Quantum 
dot 
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1. Introduction 

The Kondo effect^ in quantum dots has been an active 
research field over a decade. The early prediction about 
a characteristic gate-voltage dependence of the linear- 
response conductance^^* has been confirmed experimen- 
tally in semiconductor devices, and other new features 
of the Kondo physics are also being studied extensively in 
various situations, such as an AB ring, Joscphson junc- 
tion, ferromagnetic leads, etc. 

The equilibrium and linear-response properties of a 
single quantum dot connected to normal leads can be ex- 
plained basically based on the knowledge of the Kondo 
problem in dilute magnetic alloys,^ although there ex- 
ists some differences in experimental geometry (config- 
uration) between the impurity in magnetic alloys and 
quantum dots in semiconductor devices. Therefore, the 
low-energy properties can be described by the local Fermi 
liquid theory, ^"^^^ and the nonperturbative approaches 
developed for the alloys, such as the quantum Monte 
Carlo^'^' and numerical renormalization group (NRG) 
methods, can be applicable to the quantum dots. Par- 
ticularly, the NRG has been used successfully to calculate 
the linear-response conductance of the quantum dots.^^ 

The nonlinear transport under a finite bias voltage V , 
however, is still not fully understood, despite of a number 
of theoretical efforts. ^^"^"^ Among a variety of aspects of 
the nonequilibrium properties, in this report we focus our 
attention mainly on the low-energy properties. Specifi- 
cally, we describe the Fermi-liquid behavior of the first 
nonlinear term of the differential conductance dl/dV us- 
ing the Ward identities, which is derived by applying the 
perturbation theory in the Coulomb interaction U of Ya- 
mada and Yosida to the Keldysh Green's function. We 
show that the low-energy asymptotic form of the order 
C/^ self-energy^^ is essentially retained in all orders in U, 
and the contributions of the higher-order terms are ab- 



sorbed into the coefficients which can be written in terms 
of the local- Fermi-liquid parameters such as the width of 
the Kondo resonance A and Wilson ratio R.^^ The proof 
was provided previously in ref. 22. In the present report, 
however, we give another derivation, using the property 
of the impurity Green's function G{oj) as a functional of 
a nonequilibrium distribution function fesiijj), through 
which the dependence of G{io) on eV and T arises. This 
property also allows us to deduce some exact results in 
the limit of large eV.^^ In the present report, we re- 
examine the low-energy properties with an emphasis on 
this aspect of the Green's function as a functional. 

We also present an alternative description of the low- 
voltage Fermi-liquid behavior using the renormalized 
perturbation theory (RPT).^^ The unperturbed Green's 
function of the RPT in the Keldysh formalism consists 
of the propagators of the free quasiparticles, which are 
determined by the renormalized resonance of the width 
A and the nonequilibrium distribution /cff(w). To second 
order in the residual interaction U, which is defined in 
eq. (44), it gives the exact low-energy (eV)'^ coefficient of 
dl/dV. The higher order terms in U determine the high- 
energy properties. It has recently been confirmed that in 
equilibrium a combination of the RPT and NRG gives an 
efficient way of calculating the temperature dependence 
of the susceptibility,^^ so that the RPT seems to be one 
possible approach to the nonequilibrium properties be- 
yond the Fermi-liquid regime. 

In §2, we describe the Keldysh formalism for the An- 
derson impurity in order to describe clearly the proper- 
ties of G(w) as a functional of /off (aj). In §3, we consider 
the low-energy behavior of the self-energy at small eV 
using the Ward identities, and give an exact low-energy 
expression of dl/dV in the electron-hole symmetric case. 
In §4, the RPT is applied to the low-voltage Fermi-liquid 
regime, and the procedure of the perturbation expansion 
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in U in the Keldysh formalism is provided. In the ap- 
pendix, details of the Ward identities and properties of 
G{u)) as a functional of /cir(t^) arc given. 

2. Keldysh Formalism for the Anderson Model 

We start with the single Anderson impurity connected 
to two leads at the left (L) and right {R): 

H = He + Hd + ifmix + , (1) 



H, 



X=L,R ka 



U 



(2) 



Hd = Y^Eanda, Hu = '^[Yl^da - , (3) 
cr ^ a ' 

ifmix = E E^A (4V'A. + i'Ld. ) , (4) 

\=LM a- 

where annihilates an electron with spin a at the dot, 
nd(T = d\d^, and Ej, = Cd + U/2. We assume that the 
onsitc potential is a constant independent of the bias 
voltage, and take the Fermi level at equilibrium to 
be the origin of the energy, i.e., /it = 0. In the lead at A 
{— L, R), the energy spectrum is given by ekx = Cfe + eVA. 
To specify how the bias voltage V is applied to each of the 
leads, we introduce a parameter a\ such that = ai,V 
and Vji = —ajiV with + aj^ = 1. In cq. (4), vx is the 
tunneling matrix element between the dot and lead at A, 

and V'a.t = Sfc ^kXcrl^J^ ■ units h=\. 

In the thermal equilibrium, wc know that the density 
matrix is given by p^^ oc e~'^^, and thus the Hamil- 
tonian determines both the time evolution and statis- 
tical weight. However, in a nonequilibrium steady state 
the density matrix cannot be determined simply by if, 
and it depends on how the system has been driven to 
the steady state. The Keldysh formalism has been used 
widely for this purpose to determine the density matrix 
for nonequilibrium states. 

The method uses the procedure of an adiabatic switch- 
ing on, which is described by the operator U(t,to) = 

Texp[-i//^di'-ff2 (*')]■ Here, 0(t) = e'^^^Oe''"^* is 
an operator in the interaction representation with re- 
spect to Hi, which is a time-independent part of the 
total Hamiltonian H{t) = Hi -\- H2 e"*^!*! . In the in- 
teraction representation the density matrix defined by 
pit) = e'^i* e~'-^i* can be rewritten in the form 



p{t) = U{t, -00) p{-oo)U{-oo,t) 



(5) 



where p{—oo) represents the initial statistical weight. 
The average value of a Heisenberg operator Onit) = 
U{0, t) d{t) U{t, 0) is given by 

{Onit)) = Tr[p{0)OH{t)] 

= TV p(-oo) U{-oo, +00) U{+oo, t) d{t) U{t, -00) 

(6) 

The stream of time seen in this expression is usually il- 
histrated as the Keldysh contour shown in Fig. 1: the 
+ branch corresponds to the time evolution by the op^ 
erator ^^(-00,-^00) = Texp[i dt' i?2(i') where T 
denotes the anti-time-ordering operator. If one chooses 



—00 
—00 



-|- branch 



J 



-|-oo 



— branch 

Fig. 1. The Keldysh contour for the time evolution. 

Hi to be bilinear, the Fcynman-diagrammatic approach 
is applicable for the Green's functions defined by 

G--{t) = -i(Td,(t)4(0)), (7) 

G-+{t) = i{diiO) djt)) , (8) 

G+-(t) = -i(d,(i)4(0)), (9) 

= -i(frf,(t)4(0)). (10) 

These functions arc linearly dependent G ^ + G"' = 
G -|- Furthermore, the retarded and advanced 

Green's functions can be written as G^ = G — G 
and G" = G — G+~ , respectively. 

2.1 Traditional Formulation 

To describe a nonequilibrium steady state under a fi- 
nite bias voltage, Caroli et al?^ has introduced the initial 
statistical weight of the form 

p{-oo) oc e-^(^'^ + ^<=-''^^^-''«^«) , (11) 



The two chemical potentials. 



where = J2ka clx^CkXa 
PL = eVi and pn = cVr, are defined with respect to the 
isolated systems described by Hi = Hd + He, and the 
remaining part H2 = ifmix + Hu is switched on adiabat- 
ically. Specifically, in the noninteracting case U = 0, the 
Green's functions can be written in the form 



G^+ico) 



G+-(u;) 



[1 - .Mcu)] Gl{u) + Gliio) , (12) 

: -/eff(a;)[GS(^)-Go»(c.)], (13) 

: [l-/eff(a;)][GS(a;)-Gg((^)], (14) 

G++{w) = -[l-/eff(a;)]GS(a;) - /eff(a;) GS(a;), (15) 

where Gl{uj) = [uj- Ed + iA]"\ G'^(w) = {G'^(w)}*, 
and A = Tl -\-Tr with Fa = ttp^^a- assume that 
the density of states p^i^) = 5{uo — ^kx)/-^ ^ ^o^^' 
stant, and the band width is very large. One important 
feature we see in eqs. (12)-(15) is that all the informa- 
tion about the nonequilibrium distribution is contained 
in the distribution function, 



/eff(w) = 



(16) 



Here ^(a;) = /(w - p^), and /(w) = [e^'^ + At 
T = the distribution function fes{u)) has two steps, 
at u ~ PL and pn, as shown in Fig. 2. At eV = 0, it 
coincides with the usual Fermi function f{u)). 

The interacting Green's function G{w) satisfies the 
Dyson equation, 

{G(a;)}-i = {Go(a;)}-i-S(a;) 

Gn 



Go = 





G+- 



Gg + 

Go++ 



(17) 

(18) 
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Fig. 2. The nonequilibrium distribution fgg{ui) at T = 0. 



Here is the self-energy due to Hjj: the four elements 
are linearly dependent E + E+~ = — E — E++, 
and we have also two extra relations in the w-space, 
E»(a;) = {E''(w)}* and E" (w) = -{E++(a))}*. Using 
these relations, the retarded Green's function is written 
in the form 

Cfw) = = (19) 

with E'' = E~~ + E '". The four elements of G{ui) are 
also written in the forms similar to eqs. (12)-(15), for 
which Gq and ^re replaced by the interacting ones and 
fes{i^) is replaced by a correlated distribution defined by 



Moj) Tl + fniuj) Tfl - 1 E-+(a;) 



(20) 



This function was introduced by Hershfield et ai, and 
was studied using the order self-energy.^'' Note that 
E '"(w) is pure imaginary, and at eV = it takes the 
form E-+(w)|g^^o = 2i/(a))Im T,''{Lo)\^y^Q. Thus, in the 
nonequilibrium state, the distribution function /^(w) 
generally depends on the interaction U, while in equi- 
librium eV = it coincides with the Fermi function. 



2.2 Alternative Formulation 

As described in the above, the noninteracting Green's 
function for Hq = Hc + Hd + ffmix can be calculated ana- 
lytically taking all contributions of the tunneling matrix 
element i?mix into accoimt. A question arises: do we al- 
ways have to start with the isolated systems to obtain 
eqs. (12)-(15)? The answer is no. An alternative descrip- 
tion was given by Hershfield.'^'' The basic idea is to assign 
the two different chemical potentials directly to the left- 
and right-moving scattering states which are written for- 
mally, using the Lippmann-Schwinger equation, as 

where X = L, R. The incident wave comes in the left for 
Tfeio-' '^^ right for 7^^^. These scattering states 
are the eigenstates, by which Hq can be diagonalized as 



\=L,R ka 



(22) 



Note that generally the bound states and continuum 
states without the degeneracy (for the left and right 
movers) are present. Such states are not distinguished 
from the degenerate scattering states in eq. (22) for sim- 



plicity. With these scattering states, the density matrix 

for U = Q can be expressed explicitly as"^" 

Po(0) oc e-^( H, + H. + ff„„. - M^AT^ - mhATk ) ^ (23) 

where M\u = '^k^kXa^kXa- '^^'^ confirm that the 
noninteracting Green's functions eqs. (12)-(15) can be 
calculated directly from eq. (23). Therefore, the Coulomb 
interaction Hjj can be switched on starting from the con- 
nected system taking /5o(0), given in eq. (23), to be the 
initial statistical weight. It is carried out by using eq. 
(6) and redefining the initial condition as Hi => Hq, 
H2^ Hu, and p(-oo) ^ e''^^ " " ^ . 

The perturbation scries in U of the nonequilibrium 
Green's function G can be generated automatically using 
the path integral representation,'^^ 



(24) 



(25) 



where r]^^(t) is a Grassmann number for the branch v 
(= — , -|-) in the Keldysh contour. The action S is defined 
by 

(26) 
(27) 



So = 



/oo 
-OO 



E 



Su = 



dt 



where 

Ko{t,t/) ^ 



\lll-{t)lla-{t) 



{Go(a;)}-^-'-(*-*') 



(28) 



(29) 



and {t) = ( r]l_ {t) , ry^.,. {t) ) . In the Keldysh formalism 
the perturbation expansion works with the real frequen- 
cies (or real times). Therefore, eq. (25) shows that the 
dependence of G{ijj) on the bias voltage and tempera- 
ture arises through fes{u}) in the noninteracting Green's 
function Go{l^) which determines Sq via eqs. (27) and 
(29). Thus, the full Green's function G{co), can be re- 
garded as a functional of /eff(a;). The precise form of the 
functional is obtained by expanding e'"^ in eq. (25) in the 
power series of U, and substituting eqs. (12)-(15) into 
every single Go's in the series. Therefore, the change in 
the self-energy, SU, caused by a small variation in the 
distribution function, (5/cff, can be expressed in the form 



<5/ett=0 

higer order terms in ^i/eff 



5GZ{u') 

(30) 



where 



-[GS.,(u>')-GS.,(^')]'5./offV). (31) 
The functional derivative SHa/SGQ^, can be related to 
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the vertex corrections in the Keldysh formahsm. The 
functional aspect discussed here is analogous to the func- 
tional approach of Luttinger and Ward.^^ However, in 
eq. (30), the functional derivative is taken with respect 
to noninteracting Green's function. At finite tempera- 
tures T ^ the distribution function /eff(u;) dose not 
have the discontinuities, and thus it can be treated as 
a regular function in general discussions. Nevertheless, 
the singularities appearing in the limit of T ^ play an 
important role, for instance, as we see in eq. (39). 

3. Fermi-Liquid Behavior at Small Voltages 

In equilibrium and linear-response' rc'ginie, the low- 
energy properties at w, T <C Tk can be described by the 
local Fermi liquid theory,^" where Tk is the Kondo tem- 
perature. The Fermi liquid theory can also describe the 
nonlinear response at small bias- voltages eV -C Tk-^^ 
Our proof uses the Ward identities'^' ^^'^^ in the Keldysh 
formalism. In this section, we describe the outline of the 
derivation of the identities, and then determine the low- 
voltage behavior of the differential conductance dl/dV 
up to terms of order {eV)"^ in the electron-hole symmet- 
ric case. 

3.1 Ward identities 
We first of all consider the behavior of Go{u) at small 

is written in the form 

d 



eV. The first derivative at eV 



d{eV) 



= — a 



eV=0 



d_ 

du) dEd 



<J^O:eq(^) ) 



(32) 



where a = {aL^L — ctRTR^/iJ^L+^R)-, and the label "eq" 
in the subscript stands for the "equilibrium", so that 
<^0:e'q = G^'''|ey=o- O^ing to the properties of /eff(w), 
the differential coefficient with respect to eV can be re- 
lated to the equilibrium quantities in the right-hand side 
eq. (32). From the discussions in §2, the self-energy can 
also be regarded a functional of Go(w). Thus, the dif- 
ferential coefficients of S(a;) with respect to eV can be 
calculated taking the derivative of G?o's appearing in the 
perturbation series in U, as described in the appendix. 
Then, using eq. (32), we obtain 



d{eV) 



■ a I ^ + ^ 1 Seq(a;) , (33) 



d{eVf 



eV=0 



eV=0 



du) dEd 

d d 



+ 



duj dEd 
Tl Tr 
{TL + TRf 



Seq(a;) 



£'2Seq(a;) , (34) 



where Seq('^) = ^{^)\ev=o- '^^^ operator D"^ acts on 
the noninteracting Green's functions in the perturba- 
tion series for Seq(w), and it takes the second derivative 

{d/du' + d/dEd)'^, as 

D2Seq,^(w) = 



I ^Seq,o-(w) 



'7 + 4)'^--'(^')- 



(35) 



— 



r 



crir ik— (T 



cj (T a> cr 



— 1 








r 




cr ' 




I— cr 




r 









(T 1' ^' -o- 



co^ -cr 



to — cr 
► 



(a) (b) (c) 

Fig. 3. The diagrams contribute to the singularities. 



Using these relations, the low-bias behavior of the self- 
energy can be deduced from the equilibrium quantities. 

Specifically, at T = and eV = 0, the usual 
zero-temperature formalism is applicable for the causal 
Green's function defined with respect to the equilibrium 
ground state, 

G--{u,) = Gl^{u,)9{u,) + Gt^{u,)e{-w) , (36) 

where 0{uj) is the Heaviside step function. Thus at T = 
the causal matrix-element in the right-hand side of eq. 
(33) can be related to the vertex corrections,'^'^ 



- r,<T';<Tv(w,0;0,a;)^eq,a'(0) , (37) 

a' 

where T^rcr'-cr'ai^T^']'-^' is the vertex function for the 
causal Green's function in the T = formalism, and 
Aeq^„{u)) = — ImGgq_^(a;)/7r. Similarly, the causal ele- 
ment of £>^Seq,(7 can be written as^^ 



dto' 



{<jJ,U}';U}',Uj) Aeq,a'(w') 



w'=0 



(38) 



Eliashberg^'* has shown quite generally by using the 
Lehmann representation that the imaginary part of the 
vertex function has some singularities.'^' ^^'^^ For small 
frequencies, the singularities relevant to eqs. (37) and 
(38) arise from the diagrams shown in Fig. 3. The asymp- 
totic form of the imaginary part of eq. (38) at small u 
and u)' can be calculated as^^ 
d 



duj' 

a' 

= - |rn;iT(o,0; 0,0)1' 

duj" 



X Im 



27ri ^^'i 



-I- 



f du" 
J 2m 



G: 



(w-w' + w") 



-{U) + IV'-L0") 



■7r{Aq(o)}' |rnaT(o,0; 0,0)1' 

—2 sgn(a;' — u) + sgn(a;' -|- u) 



(39) 



Here the first term in the last line corresponds to the 

contributions of the diagram (a) and (b) , and the second 
term corresponds to that of the diagram (c). Due to these 
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singularities, the value of eq. (39) in the limit of lo, lu' 
depends on which frequency is taken first to be zero. 
Taking the limit for small frequencies carefully, the low- 
energy behavior of Im Y7'{uj) are determined up to terms 
of order cj^, {eVY , and T^, 

ImE'-(c^)= ^|{Aeq(0)}^ |rTiaT(0,0;0,0)|' 



[uj-aeVy + 



[eVf + {irTf 



(40) 



The result at equilibrium eV = has been provided by 
Yamada and Yosida,^'^ and it is extended to the nonequi- 
librium steady state here up to terms of order (eV)'^. 
Note that we have not assumed the electron-hole sym- 
metry so far. 

3.2 Results in the electron-hole symmetric case 

In this subsection we consider the low-energy behav- 
ior of G^{uj) and dl/dV using the result of lm'S^{u!) 
obtained in eq. (40). Specifically, we concentrate on the 
electron-hole symmetric case, where = — C//2, = 
Ffl, and = = 1/2. In this case Aeq(O) = l/{nA), 
and the real part of the self-energy takes the form 



ReE''(cj) = (1- 



1 - 



(41) 
(42) 



Thus, G'^Ilo) can be calculated exactly up to terms of 
order w^, and {eVY using eq. (40), 



2A(7rA)2 



(43) 



where the renormalized parameters are defined by 

A = zA, [/ = z2r||^|(0,0;0,0) . (44) 

The order result of Hershfield et al.^'^ can be repro- 
duced from eq. (43) replacing U by the bare Coulomb 
interaction U , and using the order [/^ result for the renor- 
malization factor^-'^''^^ z = 1 — (3 — 7r2/4) + • • • , where 
u = [//(ttA). 

Thus, in the symmetric case the low- voltage behavior 
is characterized by the two parameters A and U . These 
parameters are defined with respect to the equilibrium 
ground state, for which the exact Bethe ansatz results 
exist'^^"^^ as shown in Fig. 4. The width of the Kondo 
resonance A decreases with increasing U , and for u > 
2.0 it is approximated well by the asymptotic form A ~ 
{A/'k)Tk, where the Kondo temperature is defined by 



Tk = TTl\yJu/{2ii) exp[-7r2u/8 + l/(2u)] 



(45) 



The Wilson ratio is usually defined by R = Xs/l, where 
7 and Xs are the enhancement factors for the T-linear 
specific heat and spin susceptibility, respect ively.-^^ Al- 
ternatively, it can be written in terms of A and U , as^^ 

R-l ^ U/inA) . (46) 

The Wilson ratio increases with u from the noninteract- 



o ' 
tn 

■§ 0.8 

c 

o 

0.6 

oa 

■B 0.4 



o 0.2 



\ 1 1 1 


\ 


R-1 




-a/a\ V""^ 




- 








/ \-4TK/nA 






/ X 






/ 






/ 1 1 


1 — 





12 3 4 

U/7iA 

Fig. 4. The U dependence of A/A and R — I, which can also 
be interpreted as z and U /{it A), respectively. These parameters 
were calculated using the Bethe ansatz solution summarized in 
ref. 36. 



ing value i? = 1 to the strong-coupling limit value R — 2. 
The charge excitations at the impurity site are still sur- 
viving for u < 2.0, and it makes the value of R smaller 
than 2. 

The nonequilibrium current / can be calculated from 
the retarded Green's function, ■^^ 

^ rdu;[h-fn]^^^[-ImG^{io)]. 

n J_oo 1 _R + 1 L 

(47) 



Substituting eq. (43) into eq. (47), the differential con- 
ductance dl/dV can be determined exactly up to terms 
of order and [eVf, 



dl 
dV 



2e^ 
h 



1 - 



l + 2(i?- 1) 



l + 5(i?-l)2 



vrT 
X 

eV 
A 



(48) 



The result shows that the nonlinear {eVy term is also 
scaled by the resonance width A, and the coefficient gen- 
erally depends on the parameter {R~ 1)^, or U'^ /{irAy. 
As mentioned, in the strong-coupling limit u oo, the 
two characteristic parameters become A — s- (4/7r) Tk and 
R^2. 

3.3 Comparison with other approaches 

To our knowledge, similar attempts to calculate the 
coefficient cy of the {eVY term of dl/dV have been made 
by two groups^'''^^ in the strong-coupling limit u — > oo: 



dl_ 
dV 



2e^ 
h 



1 



ct 



ttT 
TV 



cv 



eV_ 
TV 



(49) 



To compare the results, the difference in a numerical 
factor of order 1 in the definition Tk must be taken 
into account. To avoid this uncertainty, we use Tk de- 
fined in eq. (45), and rescale the results presented in 
refs. 20 and 21 such that the coefficient for the linear- 
response T^ term agrees with the result of Yamada and 
Yosida ct = {tt/A)^. Kaminski, Nazarov, and Glazman^^ 
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have carried out a perturbation expansion around the 
strong- couphng fixed point to obtain = (3/8) Ct- 

Konik, Saleur, and Ludwig^^ have used the equihbrium 
Bethe ansatz solution, and then made some extra as- 
sumptions for calculating the nonlinear coefficient to ob- 
tain Cy^^ = 4:Ct (the parameters corresponding to az, 
and an used by KSL seem to be different from ours). 
Our result eq. (48), which is obtained using the Ward 
identities, shows c^'""^ = (3/2) Ct in the strong-coupling 
limit. 

Although the Hamiltonian is somewhat different, we 
also note for comparison that Schiller and Hershfield^^ 
obtained the result corresponding to Cy^ = 3ct for a 
special parameter set which can be related to the Emery- 
Kivelson solution of the two-channel Kondo model. 

4. Renormalized Perturbation Theory at Finite 
Bias Voltages 

Although the description of the low-energy properties 
discussed in §3 is exact, the underlying physics of the 
quasiparticles might not be seen directly in the micro- 
scopic derivation. In the case of the three-dimensional 
Fermi liquid, the vertex function played a central role to 
clarify a link between the intuitive picture of the quasi- 
particles and Green's functions. Specifically, the resid- 
ual interaction between two quasiparticles, which had 
been introduced phenomenologically, was shown to be 
connected to the forward scattering amplitude.'*" 

For the Anderson impurity the vertex function at 
Fermi energy, z^r||;|-f (0, 0; 0, 0), corresponds to the scat- 
ter amplitude, and it is equal to U by the definition in eq. 
(44). The perturbation expansion in U, which has been 
formulated precisely in the equilibrium case by Hew- 
son,^^ provides the link between the quasiparticles and 
microscopic theory of the local Fermi liquid. All the basic 
Fermi- liquid behavior have been shown to be reproduced 
in the expansion up to terms of order U^. Furthermore, 
the approach is not limited to low energies. To carry out 
the expansion systematically, however, one has to take 
account of the renormalization conditions that are nec- 
essary to avoid overcounting of the many-body effects, 
because in the renormalized perturbation theory (RPT) 
the expansion parameter already contains some contri- 
butions of the Coulomb interaction. 

In this section we apply the RPT to the nonequilib- 
rium steady state. It reproduces the result of dl/dV in 
the Fermi-liquid regime, and gives us one possible way 
to calculate the corrections needed at high voltages. For 
simplicity, we concentrate on the electron-hole symmet- 
ric case; = —U/2, Tl = Tr, and oil = otR = 1/2. 
The unperturbed Green's function is defined such that 
it describes the Kondo resonance with the renormal- 
ized level width Gq{uj) = -|- iA] , as that in the 
equilibrium case.^^ However, in the nonequilibrium case, 
it is not obvious how the distribution function for the 
free quasiparticles should be given by. We simply as- 
sume here that it is given by the noninteracting one, 
which in the electron-hole symmetric case takes the form 
f^g{uj) = [/(w - eV/2) + f{uj + eV/2) ] /2. Hence, the 
four elements of the unperturbed Green's functions. Go, 




Fig. 5. The second-order diagram of S. The dashed represents U, 
and solid lined represents the free-quasipaxticle propagator Go- 



take the forms 
Go~(a;) = [1 - fM] GU^) + fes{co) Gg(a;), (50) 



G^+{u;)^-fM GS(c.)-Gg(c.) 



G+-{u;) = [1 - /eff(a;)] GS(u;) - GS(a;) 



(51) 
(52) 



Go++(a;) = - [1 - /eff(a;)] Gg(w) - /eff(w) G5(a;) . (53) 

Correspondingly, the full propagator of the quasiparti- 
cles, which includes all effects of U, is defined by G{co) = 
z~^G{lj). Therefore, in terms of the renormalized quan- 
tities, the nonequilibrium current, eq. (47), is written as 



_ 2e r 



(54) 



The self-energy correction due to U satisfies the Dyson 
equation of the form S(w) = {Go(w)}~* — {G{uj)}~^. 
Using eq. (17), we have 

S(a;) = {Go(a;)}-i - z [{Go{ij)}-' - 5](a;)] 

= zi:{uj) + {l-z)T3iO , (55) 
where Tj for i = 1, 2, 3 is the Pauli matrix. 

4-1 Low-energy behavior up to terms of order u^, 
and (el/)2 

At low-energies, the ti;-linear contributions in the right- 
hand side of eq. (55) cancel out owing to eq. (41). The 
contributions of order ui^, and (eV)^ arise from the 
second-order diagram for S(a;) shown in Fig. 5, where 
the solid and dotted lines represent Go{lo) and U, respec- 
tively. Calculating the contributions from the diagram, 
and then taking the cancellation of the w-linear term into 
account, we have 



S'-(a;) 



2A{wAf 



I [eVf + {nTf 



+ 



(56) 

It simply reproduces the renormalized Green's function 
G''(w) corresponding to eq. (43). Furthermore, using eq. 
(54), the (eV)^ term of dl/dV in eq. (48) is also repro- 
duced in the expansion up to terms of order U"^. Since 
eqs. (43) is asymptotically exact, the higher-order terms 
in U do not change the low-energy behavior in eq. (56). 

Note that eq. (56) follows from the fact that we have 
used /eff(a;) for the distribution function of the free quasi- 
particles. This assumption seems to be justified also from 
the fact that the many-body effects on the correlated dis- 
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tribution f^ffito) defined in eq. (20) arise first in the order 
contribution. 

4..2 Beyond the Fermi-liquid regime 

To study the higher-energy behavior at large lo, T, and 
eV with the RPT, one needs to calculate the higher order 
terms in U. In the following, we describe the outline of 
the procedure of the expansion. 

At high energies the renormalization factor z cannot 
be defined with respect to T = and eV = no longer. 
This is because the coefficient of the a)-linear term of 
the self-energy depends on T and eV . For instance, in 
the next order, the terms of the form T^w and {eVY'uj 
exist. Therefore, z is redefined such that the w-linear 
contributions in eq. (55) cancel out 



as''(w) 



dw 



= 



(57) 



Hence z generally depends on T and eV. The perturba- 
tion expansion in U can be carried out following that in 
the equilibrium case.^^ We first of all rewrite the action 
S in the form 

S = z-^So + z-^SuU/U - 5cou (58) 
Scon = A z-^SuU/ U + z-^So - So , (59) 

where Sq is the action for the free quasiparticlc corre- 
sponding to the propagator Go(a;), and A = 1 — z"^ U/U. 
In eqs. (58) and (59) the factor is introduced just 
to cancel the bare Coulomb interaction U included in 
Su by the definition in eq. (28). The perturbation se- 
ries in U is generated by taking z~^So in eq. (58) to 
be the unperturbed part and taking the remaining terms 
z~'^SuU/U — Scou to be the perturbed part. Here, 5cou is 
the counter-term which avoids overcounting of the many- 
body effects. Specifically, the last two terms in the right- 
hand side of eq. (59), which can be rewritten in the form 

/oo 

(60) 

It corresponds to the counter-term for the renormaliza- 
tion factor z. In the RPT, the two parameters z and A 

are regarded as functions of the rcnormalizcd parameters 
A and U, and are expanded as series in the powers of U. 
Then the expansion coefficients for z and A can be de- 
termined such that the two renormalization conditions, 
eqs. (57) and (61), are satisfied by each order in U] 



ruaT(o,0;0,o) = u . 



(61) 



Here F 



T-l;iT 



z^F|m| is the vertex part for the four 
external causal Green's functions G~~ . In the RPT, 
I^Ti;iT(0;0;0, 0) is calculated in the power series in [/, 
and at high-energies it generally depends on T and eV . 
Note that the^ontribution of the parameter A first arises 
in the order terms. For this reason, the condition 
of A is not necessary to be taken into account in the 
expansion up to order U^. _ 

As already mentioned, higher-order terms^in U are 
needed to study the high-energy behavior of S'"(a;) and 



dl/dV beyond the u)^- and (eF)^-terms. One possibility 

is to include the contributions up to terms of order f/^. 
The corresponding calculations in the bare-?7 expansion 
have been carried out by Fujii and Ueda.^^ Alternatively, 
in the equilibrium case at eV = 0, a combination of the 
RPT and NRG has been examined recently, and the re- 
sults reproduce the T-dependence of the spin suscepti- 
bility accurately in a wide temperature range. Such a 
combination would be another possibility to go beyond 
the Fermi-liquid regime at large bias voltages. 

5. Summary 

We have studied the low-energy properties of the An- 
derson model under a finite bias voltage V using the 
properties of the Kcldysh Green's function at the impu- 
rity site G{uj) as a functional of the nonequilibrium dis- 
tribution function /cii (cl)). Through the distribution func- 
tion /eff(w), the T- and eF-dependence of G(a;) arise. 
The Ward identities for the derivative of the self-energy 
with respect to eV follow from these properties that can 
be summarized in the form of eq. (30). Using the Ward 
identities, the differential conductance dl/dV has been 
determined up to terms of order (eV)'^ in eq. (48) in the 
electron-hole symmetric case. The coefficients are deter- 
mined by two characteristic parameters A and R. 

We have also described the low-energy properties us- 
ing the renormalized perturbation theory in the Keldysh 
formalism. To second order in U, it reproduces the exact 
(eV)"^ coefficients for dl/dV. The Fermi-liquid behavior 
of dl/dV follows from the assumption that the distribu- 
tion func;tion for the free quasiparticles are the same as 
that of the noninteracting electrons /efi(w). In order to 
study the corrections to the Fermi liquid theory at large 
bias voltages with the J^PT, one needs to calculate the 
higher order terms in U . 
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Appendix: Derivation of eqs. (33)— (34) 

The first derivative of the self-energy with respect to 
eV can be written, using eq. (30), as 

eV=0 



d{eV) 



doj' 



(AT) 



eV=0 



Similarly, the derivative of ^^(a') with respect to u) at 
eV — is written in the form, 

d_^_d_ 

diO 



BE, 



Seq,<7('^) 



dw' 



(5Seq,a(w) 



d d 



dw dE^ 



GZ^,a'{^' + ^) 
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d 

duj' 



dE„ 



^0:eq,cr 



\J). (A.2) 



Here we have used the property that the frequency lo' can 
be shifted to w' + a; without changing the result. This is 
because the vahic of the self-energy docs not change if all 
the frequencies which are assigned to the Green's func- 
tions in a closed-loop diagram are shifted by the same 
amount. Therefore, substituting eq. (32) into eq. (A-1), 
we obtain eq. (33). 

To calculate the second derivative, the variation of the 
self-energy (5Scr must be calculated up to terms of order 
(<5/cfl)^, and we find 



52S,(w) 



d{eVf 

= E 



ey=o 



dw' 



eV=0 



dw'do; 



<5'Seq,a(w) 



d{eY) d{eV) 



eV=0 

(A-3) 



Then flowing along the similar line, we obtain eq. (34) 
using cq. (32) and the corresponding relation for the sec- 
ond derivative 

eV=0 



d{evy 



where k = (air^ + al,rR)/{TL + Tr). 
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